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Abstract 

In a previous work [E.M. Prodanov, R.I. Ivanov, and V.G. Gueorguiev, Reissner- 
Nordstrom Expansion, Astroparticle Physics 27 (150-154) 2007], we proposed a 
classical model for the expansion of the Universe during the radiation-dominated 
epoch based on the gravitational repulsion of the Reissner-Nordstrom geometry — 
naked singularity description of particles that "grow" with the drop of the tempera- 
ture. In this work we model the Universe during the Reissner-Nordstrom expansion 
as a van der Waals gas and determine the equation of state. 
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1 Introduction 



In 1971, Hawking suggested [T] that there may be a very large number of gravitation- 
ally collapsed charged objects of very low masses, formed as a result of fluctuations in 
the early Universe. A mass of 10 14 kg of these objects could be accumulated at the 
centre of a star like the Sun. Hawking treats these objects classically and his arguments 
for doing so are as follows pQ: gravitational collapse is a classical process and micro- 
scopic black holes can form when their Schwarzschild radius is greater than the Planck 
length (G/i/c 3 ) -1 / 2 ~ 10~ 35 m (at Planck lengths quantum gravitational effects do not 
permit purely classical treatment). This allows the existence of collapsed objects of 
masses from 10 -8 kg and above and charges up to ±30 electron units p]. Additionally, 
a sufficient concentration of electromagnetic radiation causes a gravitational collapse 
- even though the Schwarzschild radius of the formed black hole is smaller than the 
photon's Compton wavelength which is infinite. Therefore, when elementary particles 
collapse to form a black hole, it is not the rest Compton wavelength hc/mc 2 that is to 
be considered — one should instead consider the modified Compton wavelength hc/E, 
where E ~ kT >>> mc 2 is the typical energy of an ultra-relativistic particle that went 
to form the black hole pQ. Microscopic black holes with Schwarzschild radius greater 
than the modified Compton wavelength hc/E, can form classically and independently 
on competing quantum processes. 

Hawking suggests that these charged collapsed objects may have velocities in the range 
50 - 10000 km/s and would behave in many respects like ordinary atomic nuclei pQ. 
When these objects travel through matter, they induce ionization and excitation and 
would produce bubble chamber tracks similar to those of atomic nuclei with the same 
charge. The charged collapsed objects survive annihilation and, at low velocities (less 
than few thousand km/s), they may form electronic or protonic atoms pQ: the positively 
charged collapsed objects would capture electrons and thus mimic super-heavy isotopes 
of known chemical elements, while negatively charged collapsed objects would capture 
protons and disguise themselves as the missing zeroth entry in the Mendeleev table. 
Such ultra-heavy charged massive particles (CHAMPS) were also studied by de Rujula, 
Glashow and Sarid [2] and considered as dark matter candidates. 

Dark Electric Matter Objects (DAEMONS) of masses just above 10~ 8 kg and charges of 
around ±10 electron units have been studied in the Ioffe Institute and positive results 
in their detection have been reported [3] — observations of scintillations in ZnS(Ag) 
which are excited by electrons and nucleons ejected as the relic elementary Planckian 
daemon captures a nucleus of Zn (or S). 

The DAMA (DArk MAtter) collaboration also report positive results [3] in the detection 
of such particles using 100 kg of highly radiopure Nal(Tl) detector. 
Such heavy charged particles can serve as driving force for the expansion of the Uni- 
verse during the radiation-dominated epoch in a classical particle-scale model, which 
we recently proposed |5 . Along with this type of particles, within our model, magnetic 
monopoles can also play the same role for the expansion of the Universe: it has been 
suggested [6] that ultra-heavy magnetic monopoles were created so copiously in the 
early Universe that they outweighed everything else in the Universe by a factor of 10 12 . 
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Our particle-scale model gives the expected prediction for the behaviour of the scale fac- 
tor of the radiation-dominated expanding Universe, a(r) ~ yfr , and can be considered 
as a complement to the large-scale Friedmann-Lemaitre-Robertson-Walker (FLRW) 
model (see, for example, [TJ [8]) which describes the Universe as isotropic and homoge- 
neous, with very smoothly distributed energy-momentum sources modeled as a perfect 
fluid, applicable on scales much larger than galactic ones. 

This recently proposed [5 J classical mechanism for the cosmic expansion models the Uni- 
verse as a two-component gas. One of the fractions is that of ultra- relativistic "normal" 
particles of typical mass m and charge q with equation of state of an ideal quantum 
gas of massless particles. The other component is "unusual" — these are the particles 
of ultra-high masses M (of around 10 -8 kg and above) and charges Q (of around ±10 
electron charges and above) — exactly as those described earlier. 

For an elementary particle such as the electron, the charge-to-mass ratio is q/m ~ 10 21 
(in geometrized units c = 1 = G), while for the "unusual" particles, M<Q. In view 
of this, the general-relativistic treatment of elementary particles or charged collapsed 
objects of very low masses also necessitates consideration from Reissner-Nordstrom (or 
Kerr-Newman) viewpoint — for as long as their charge-to-mass ratio remains above 
unity. We also treat the "unusual" particles classically (in line with Hawking's argu- 
ments outlined earlier). That is, the "unusual" particles are modelled as Reissner- 
Nordstrom naked singularities and the expansion mechanism is based on their gravito- 
electric repulsion. Instead of the Schwarzschild radius, the characteristic length that is 
to be considered now and compared to the modified Compton length [I], will be the 
radius of the van der Waals-like impenetrable sphere that surrounds a naked singularity 
(sec |9j for a very thorough analysis of the radial motion of test particles in a Reissner- 
Nordstrom field). As shown in [S], for temperatures below 10 31 K, the radius of the 
impenetrable sphere of an "unusual" particle of mass 10 -8 kg and charge ±10 electron 
units is greater than the modified Compton wavelength of the " unusual" particle itself. 
Naked singularities have been subject of significant scrutiny for decades. In the 1950s, 
the Reissner-Weyl repulsive solution served as an effective model for the electron. Very 
recently, a general-relativistic model for the classical electron — a point charge with 
finite electromagnetic self-energy, described as Reissner-Nordstrom (spin 0) or Kerr- 
Newman (spin 1/2) solution of the Einstein-Maxwell equations, — has been studied by 
Blinder [10]. Naked singularities are disliked — hence the Cosmic Censorship Conjecture 
|11] — but not ruled out — there is no mathematical proof whatsoever of the Cosmic 
Censorship. At least one naked singularity is agreed to have existed — the Big Bang 
- the Universe itself. Of particular importance in the study of naked singularities are 
the work of Choptuik [12] , where numerical analysis of Einstein-Klein-Gordon solutions 
shows the circumstances under which naked singularities are produced, and the work 
of Christodoulo [13] who proved that there exist choices of asymptotically flat initial 
data which evolve to solutions with a naked singularity. The possibility of observing 
naked singularities at the LHC has been studied in p3] — for example, a proton-proton 
collision could result in a naked singularity and a set of particles with vanishing total 
charge or with one net positive charge — an event probably undistinguishable from 
ordinary particle production. In a cosmological setting, naked singularities have been 
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well studied and classified — see, for example, [F5\ . 



2 Reissner— Nordstrom Expansion 

Consider a "normal" particle of specific charge q/m, and an "unusual" particle of charge 
Q such that sign(Q)q/m > — 1. If the "normal" particle approaches the "unusual" 
particle from infinity, the field of the naked singularity is characterized by three regions 
0E]: 



U(r) 




Figure 1: The potential of the interaction between an "unusual" particle of mass 
M and charge Q > and a test particle of specific charge q/m > — 1. If the specific 
charge of the probe is smaller than —1, then, as shown in 9 ], the probe will reach 
the singularity. Note that electrically neutral probes, in addition to the attraction, 
also suffer repulsion, while probes of specific charge q/m > M/Q are always 
repelled (the gravitational attraction cannot overcome the electric repulsion). The 
form of this potential is derived later Q in this section. The vertical asymptote 
to each graph is at r = ro (T) — the radius of the van der Waals-like impenetrable 
sphere surrounding the naked singularity. When minima are present, they are 
located at the critical radius r — r c — where attraction and repulsion interchange 
— see equation Q. 

(a) Impenetrable region — between r = and r = tq(T). 

For an incoming test particle, the condition for reality of the kinetic energy leads 
to the existence of two turning radii [5, 9 with a forbidden region in-between. The 
upper (outside) radius, which we denote ro(T), can be thought of as a radius of 
an "impenetrable" sphere surrounding the naked singularity. It depends on the 
energy of the incoming particle (or the temperature T of the "normal" fraction 
of the Universe): the higher the energy (or the temperature), the deeper the 
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incoming particle will penetrate into the gravitationally repulsive field of the naked 
singularity. 

(b) Repulsive region — between the turning radius tq{T) and the critical radius 
T c > r {T). 

The critical radius r c is where the repulsion and attraction interchange (we deter- 
mine r c later in this section). As the temperature drops, the "unusual" particles 
"grow" (incoming particles have lower and lower energies and turn back farther 
and farther from the naked singularity). When the temperature gets sufficiently 
low, the radius of the "unusual" particles Tq(T) grows to r c (but not beyond r c , as 
the region r > r c is characterized by attraction and an incoming particle cannot 
turn back while attracted). This means that incoming particles have such low 
energies that they turn back immediately after they encounter the gravitational 
repulsion. Incoming particles of charge q such that qQ > Mm do not even experi- 
ence attraction — we shall see that the repulsive region for such particles extends 
to infinity (the gravitational attraction will not be sufficiently strong to overcome 
the electrical repulsion). 

(c) Attractive region — from the critical radius r c to infinity. Again, there is no 
gravitationally attractive region for an incoming particle such that qQ > Mm. 

As shown in [9 , when an incoming particle has sufficiently large charge which is also 
opposite in sign to that of the naked singularity: sign(Q)q/m < —1 , the particle will 
collide with the naked singularity. When the naked singularity "captures" such parti- 
cle, its charge Q decreases and its mass M increases. If sufficient number of incoming 
particles are captured, Q will eventually become equal to M — the naked singularity 
will pick a horizon and turn into a black hole. This black hole will evaporate quickly af- 
terwards. We will assume that our "unusual" particles have survived such annihilation. 
We will also assume that these super-heavy charged particles have survived annihilation 
through all other different competing mechanisms — for example, they could recombine 
into neutral particles or decay before or after that (see Ellis et al. [16 on the astrophys- 
ical constraints on massive unstable neutral relic particles and Gondolo et al. [IT] on 
the constraints of the relic abundance of a dark matter candidate — a generic particle 
of mass in the range of 1 — 10 14 TeV, lifetime greater than 10 14 — 10 18 years, decaying 
into neutrinos). 

An interesting general-relativistic effect (with no classical analogue) is related to the 
ability of naked singularities to capture probes of charge having the same sign. This is 
associated with the inner turning radius which we denote by po{T). On Figure 2, the 
curves representing the two turning radii, ro(T) and po(T), are given as functions of 
the specific charge g/m of the probe for different temperatures. The forbidden region 
is between the two curves. As can be seen, the lower curve po{T) corresponds to a 
turning radius (capturing) of a radially outgoing probe with charge having the same 
sign as the centre. This has no classical analogue and we argue that it could serve as a 
possible mechanism for the formation of the " unusual" particles in the extremely dense 
very early Universe. Moreover, this can also allow the extension of the range of validity 
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of our model to account for the inflation of the Universe: if charge non-conservation of 
the naked singularities occurs (naked singularities picking up charge), then accelerated 
expansion can be achieved: a(r) ~ e Hr or a(r) ~ r n , with n > 1. 




Figure 2: Plotted as functions of the specific charge q/m of a radially moving test 
particle in the field of a positively charged naked singularity, the two temperature- 
dependent curves tq{T) and po(T) represent the outer and inner turning radii, 
respectively. The region between the curves ro(T) and po(T), for any given tem- 
perature T, is not allowed as it is characterized by negative kinetic energy. The 
explicit form of these curves is given in |5| later in this section (see also [HE]). 

Our expansion model assumes that initially, at extremely high energies and pressures of 
the very early Universe, the "normal" particles are within the gravitationally repulsive 
regions of the "unusual" particles with radial coordinates just above the upper turning 
radius tq{T). The particles from the "normal" fraction "roll down" the gravitation- 
ally repulsive potentials of the "unusual" particles and in result the Universe expands. 
The addition of a new class of particles (the " unusual" ) in the picture of the Universe 
does not challenge our current understanding of the physical laws governing the Uni- 
verse. The "unusual" particles interact purely classically with the "normal" component 
of the Universe and this classical interaction results in the appearance of a repulsive 
force. Our aim is to offer a possible explanation for the expansion of the Universe 
while conforming with the well established theoretical models. As shown in [5], during 
the Reissner-Nordstrom expansion, the standard relation between the scale factor of 
the Universe a and the temperature T holds: aT = const. Also, during the Reissner- 
Nordstrom expansion, the time-dependance of the scale factor is: a(r) ~ yfr (see [5] for 
details). Such is the behaviour of the scale factor during the expansion of the Universe 
throughout the radiation-dominated era, obtained by the standard cosmological treat- 
ment. 

On a large scale, the Universe is isotropic and homogeneous and for a FLRW Universe 
(0 |8]), the energy-momentum sources are modeled as a perfect fluid, specified by an 
energy density and isotropic pressure in its rest frame. This applies for matter known 
observationally to be very smoothly distributed. On smaller scales, such as stars or 
even galaxies, this is a poor description. In our picture, the Universe has global FLRW 
geometry, but locally it has Reissner-Nordstrom geometry. The compatibility of local 
Reissner-Nordstrom geometry with global FLRW geometry has been well established: 
in 1933, McVittie [18] proposed a metric embedding a Schwarzschild solution [19 in a 
FLRW universe. In 1993, Kastor and Traschen (KT) [20] found a solution desribing a 



6 



system of an arbitrary number of charged black holes in the background of a de Sitter 
universe |21| . The case of vanishing cosmological constant in the KT solution corre- 
sponds to the static Majumdar-Papapetrou (MP) solution [22], while the solution with 
positive cosmological constant is highly dynamical and describes black holes exchanging 
radiation with the background until becoming extreme (\Q\ = M). A spinning version 
of the MP solution with naked singularities was found in [23] and [23]. In 1999, the 
KT solution was extended |25j to multi-Kerr-Newman-de Sitter black holes. Metric for 
Reissner-Nordstrom black holes in an expanding/contracting FLRW universe was ob- 
tained in [26 . The interplay between cosmological expansion and local attraction in a 
gravitationally bound system is studied in [27] where new exact solutions are presented 
which describe black holes perfectly comoving with a generic FLRW universe. 
Returning to the local Reissner-Nordstrom geometry, on the level of the interaction be- 
tween the "unusual" particles and the "normal" particles of the Universe, the density and 
pressure variables should be different from those used in the large-scale geometry. We 
are going to complement the entire radiation-dominated era with Reissner-Nordstrom 
expansion and model the interaction between the "unusual" particles and the "normal" 
particles as interaction between the components of a van der Waals gas. Modeling the 
Universe as a van der Waals phase is possible in the light of the deep analogies between 
the physical picture behind the Reissner-Nordstrom expansion and the classical van 
der Waals molecular model: atoms are surrounded by imaginary hard spheres and the 
molecular interaction is strongly repulsive in close proximity, mildly attractive at inter- 
mediate range, and negligible at longer distances. The laws of ideal gas must then be 
corrected to accommodate for such interaction: the pressure should increase due to the 
additional repulsion and the available volume should decrease as atoms are no longer 
entities with zero own volumes (see, for example, [28 J. 

As an interesting development in a similar vein, one should point out the work [29 (sec 
also the references therein) which studies van der Waals quintessence by considering 
a cosmological model comprising of two fluids: baryons, modelled as dust (large-scale 
structure fluid) and dark matter with a van der Waals equation of state (background 
fluid). Van der Waals equation of state for ultra- relativistic matter has been studied by 
[30]. 

During the Reissner-Nordstrom expansion, once the temperature drops sufficiently low 
so that ro(T) becomes equal to r c , the "normal" particle with charge q, such that 
s\gn(Q)q/m > —1 and also qQ < mM, will be expelled beyond r = r c (as vq(T) < r 
always) — into the region of gravitational attraction. Due to its ultra-high energy, the 
"normal" particle will overcome the gravitational attraction and will escape unopposed 
to infinity. Thus the gravitationally attractive region is of no importance for such par- 
ticles and for them we can assume that the potential of the naked singularity is infinity 
from r = to r = ro(T) and zero from r = r c to infinity. 

For "normal" particles such that qQ > mM, the potential gradually drops to zero to- 
wards infinity (there is no attraction for these probes). For ultra-high temperatures, the 
energy E of a "normal" particle is of the order of kT. At temperatures below 10 10 K, the 
dominant term in the energy E becomes the particle's rest energy mc 2 (throughout the 
paper we use geometrized units) and, as we shall see, the turning radius ro(T) becomes 
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infinitely large below such temperature. As we model the entire radiation-dominated 
epoch with Reissner-Nordstrom repulsion, at Recombination (the end of this epoch: 
trecomb ~ 300 000 years), the free ions and electrons combine to form neutral atoms 
(q = 0) and this naturally ends the Reissner-Nordstrom expansion — a neutral "nor- 
mal" particle will now be too far from an "unusual" particle to feel the gravitational 
repulsion (the density of the Universe will be sufficiently low). During the expansion, 
the volume V of the Universe is proportional to the number N of "unusual" particles 
times their volume (one can view the impenetrable spheres of the naked singularities 
as densely packed spheres filling the entire Universe). At Recombination, V ~ trecomb ■ 
Therefore, at Recombination, the radius ro(T) of an "unusual" particle will be of the 
order of R c = N~ 1 ^t recom b . During the expansion, a "normal" particle is never farther 
than Tq(T) from an "unusual" particle. We will request that once Tq(T) becomes equal 
to R c = N~ l /^t recom } } , then the potential of the interaction between a naked singularity 
and a particle of charge q, such that qQ > mM, becomes zero. 

In this paper we use a standard treatment [28] to model the van der Waals phase of the 
Universe as a real gas and, using the virial expansion, we obtain the gas parameters. 
Combining the van der Waals equation with aT = const, we find the equation of state 
describing the classical interaction between the ordinary particles in the Universe and 
the "unusual" particles. 

Consider the Reissner-Nordstrom geometry |31[ [8] in Boyer-Lindquist coordinates {32] : 

ds 2 = - ^ dt 2 + r —dr 2 + r 2 d6 2 + r 2 sm 2 9 d<p 2 . (1) 
r A 

where: A = r 2 — 2Mr + Q 2 , M is the mass of the centre, and Q — the charge of the 
centre. We will be interested in the case of a naked singularity only, namely: Q > M. 
The radial motion of a test particle of mass m and charge q in Reissner-Nordstrom 
geometry can be modeled by an effective one-dimensional motion of a particle in non- 
relativistic mechanics with the following equation of motion [9] (see also [33] for 
Schwarzschild geometry) : 



2 + 



x _s_q e) M i, 1 _s\9t] = £^l ( 2 ) 

m M ) r 2 V m 2 ) r 2 1 " 2 ' 1 ' 



where e = E/m is the specific energy (energy per unit mass) of the three-dimensional 
motion. The expression in the square brackets is the effective non-relativistic one- 
dimensional potential and the specific energy of the effective one-dimensional motion is 
(1/2) (e 2 — 1). As we will not be interested in the effective one-dimensional motion, we 
will proceed from equation ^ to derive an expression that will serve as gravitational 
potential energy U (r) of the three-dimensional motion. In the rest frame of the probe 
(r = 0), equation ^ becomes a quadratic equation for the energy e. The bigger root 
of this equation is exactly the gravitational potential energy U (r) plus the rest energy 
m (see also [S])- Namely: 



qQ + mVA qQ 2M Q 2 
U{r) = m = — — h my 1 — | ^ — m ■ (3) 
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Since M ~ Q ~ 10 34 cm, expression ([3]) for the potential energy U(r), for distances 
above 10 -34 cm, can be approximated by: 



TT . . mM qQ m, , ~ _ 9s 1 

^a ^ ^a ^_ 



(4) 



From now on, we will use this pseudo-Newtonian potential to mimic general-relativistic 
effects with a classical theory. 

Motion is allowed only when the kinetic energy is real. Equation Q determines the 
region (r_ , r + ) within which motion is impossible. The turning radii are given by [51 E]: 



M 



1 



m M 



q_Q_ 
m M 



L) 2 -(l- e 2)(] 



m 



2 J M 2 



(5) 



We identify the impenetrable radius vq(T) of an "unusual" particle as the bigger root r + 
and the inner turning radius po(T) as the smaller root r_. The expansion mechanism is 
based on the fact that Vq(T) is inversely proportional to the temperature, namely, the 
naked singularity drives apart all neutral particles and particles of specific charge q/m 
such that sign(Q)q/m > — 1. 

Note that when e — * 1 (which happens when the rest energy becomes the dominant 
term, i.e. when kT drops below m, or below 10 10 K), then the turning radius ro(T) 
tends to infinity. 

At the point where gravitational attraction and repulsion interchange, there will be no 
force acting on the incoming particle. That is, this is the point where the derivative of 
the potential Q vanishes: 



M 



\M 2 J\ m.MJ 



(6) 



Obviously, the critical radius r c for an incoming particle charged oppositely to the 
"unusual" particle (qQ < 0) will be smaller than the critical radius for a neutral (q = 
0) incoming particle (neutral particles suffer repulsion) as the region of gravitational 
repulsion will be reduced by the additional electrical attraction. When the incoming 
probe has charge with the same sign as that of the "unusual" particle and qQ > mM, 
then r c does not exist. This means that there will be a region of repulsion only - 
the gravitational attraction will not be sufficiently strong to overcome the electrical 
repulsion. 

Finally, the potential energy of a charged probe in the field of an "unusual" particle can 
be written as follows: 



U(r) 



oo , 



r < r (T) 



-^ + f + f(-M 2 + Q 2 ) 4 , r (T)<r<R , 



(7) 



I o, 



r > R, 
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where: 

J r c , sign(Q)q/m > -1 and qQ < mM , . 
\ Rc, qQ> mM. 1 j 

Obviously, the expansion beyond r c will be due to those particles that satisfy qQ > mM. 



3 Van der Waals Equation of State 

Next, we consider the thermodynamics of a real gas. The virial expansion relates the 
pressure p to the particle number N, the temperature T and the volume V [28] : 



NkTr 



P 



V 



1 + —F(T) + (-) G(T) + 



N\2 



0) 



where the correction term F(T) is due to two-particle interactions, the correction term 
G{T) is due to three-particle interactions and so forth. We will ignore all interactions 
involving more than two particles. The correction term F{T) is |28j: 



F(T) = 2vr / A(r) r 2 dr = /3 



a 

kf ' 



(10) 



where A(r) is given by: 

A(r) = 1 

Then "van der Waals" equations is |28j : 



P+ ( 



N\2 



V 



a 



NkT / N n 
■(1 + T7/3 



V 



V 



(11) 



(12) 



In the limit N(3/V — * 0, this equation reduces to the usual van der Waals equation |28j : 



/N\2 i i 

P+{y) «](l 



N 
V 



P 



NkT 
V 



(13) 



We now assume that the "unusual" particles leave "voids" in the Universe where "nor- 
mal" particles cannot enter. Thus, the effective space left for the motion of the "normal" 
component of the gas is reduced by N/3, where (3 is the "volume" of an "unusual" par- 
ticle and N is the number of "unusual" particles. We will also pretend that "unusual" 
particles are not present and that the potential in which the "normal" particles move is 
not due to the "unusual" particles, but rather to the two-particle interactions between 
the "normal" component of the gas. In essence, we "remove" N "unusual" particles out 
of all particles and we are dealing with a gas of n "normal" particles. The "van der 
Waals" equation (12) then becomes: 



N\2 



p + (v) 



nkT , N ^ 
a=—(l + v (3), 



(14) 
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For the potential determined in ([7]), we have: 

r i. 



A(r) 



V(r) 

1 — e kr 



U{r) 
kT ' 



0, 



r < r (T) , 

ro(T) <r<R, 
r > R. 



(15) 



We then get: 



a 



r (T) 

2ir / r 2 dr 



2tt 



o 



3 r 3 (T) = -v (T) 



2tt J U{r)r 2 dr = vrmM 2 (l - ^ ) [fl - r (T)] 



(16) 



ro(T) 



irmM 



(17) 



where vq(T) is the "volume" of an "unusual" particle. Note that both a and (3 depend 
on the temperature via the particle's radius tq(T). 

We have shown [5] that for our expansion model, the standard relation between the 
scale factor of the Universe a and the temperature T holds: aT = const. Let p denote 
the density of the Universe. Then, as the volume V of the Universe is proportional to 
the third power of a and as U ~ 1/p, we have T ~ p 1 / 3 . Therefore, T/V ~ /> 4//3 . 
The volume V of the Universe during the van der Waals phase is proportional to the 
volume vq{T) of the "unusual" particles times their number N. Using equation (16), 
namely: /3 = ^vq(T), it immediately follows that N/3/V is, essentially, constant. 



Equation ( 14 ) is the equation of state for the van der Waals phase of the expanding 



Universe and can be written: as: 



P = W 4/3 " ^ 



(18) 



Here r\ is some constant. The second term depends on the temperature via a and j3 
and becomes irrelevant towards the end, as a — > when vq{T) — > R. Note also that the 
correction term —a/f3 2 is positive as a is negative. 
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